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Abstract 

Using the concept of finite-size scaling, Monte Carlo calculations of various models 
have become a very useful tool for the study of critical phenomena, with the system 
linear dimension as a variable. As an example, several recent studies of Ising models 
are discussed, as well as the extension to models of polymer mixtures and solutions. 
It is shown that using appropriate cluster algorithms, even the scaling functions 
describing the crossover from the Ising universality class to the mean-field behavior 
with increasing interaction range can be described. Additionally, the issue of finite- 
size scaling in Ising models above the marginal dimension {d* = 4) is discussed. 
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1 Introduction 

It is a common belief that at the present time, about 30 years after the 
renormalizat ion-group theory of critical phenomena was invented [1], static 
critical behavior of systems in thermal equilibrium is rather well understood. 
In particular, this is expected to be true for the most intensively studied 
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case, the Ising universality class [2], to which systems such as uniaxial ferro- 
magncts, binary alloys, simple fluids, fluid mixtures, polymer solutions and 
polymer blends belong [3]. However, in the present work, we shall draw at- 
tention to some aspects of critical behavior in Ising-like spin systems which 
are, even today, still incompletely understood. The first of these concerns the 
problem of crossover between the Ising universality class and mean-field crit- 
ical behavior. This crossover occurs, for instance, when the interaction range 
(and hence the "Ginzburg number" G entering the Ginzburg criterion [4]) is 
varied [5-11]. A closely-related crossover is found for symmetrical polymer 
mixtures when the chain length N of the polymers is varied [3,12-21]. A part 
of this crossover (though typically not the full extent of the crossover scal- 
ing function) can be probed experimentally near the critical point of fluids 
and fluid binary mixtures [22-24]. While the Ginzburg criteria [4,13,14] pro- 
vide a qualitative understanding of this crossover, the quantitatively accurate 
theoretical prediction of the crossover scaling function is a challenging prob- 
lem [25-33], and hence Monte Carlo studies [5-11,15-17] are of great potential 
benefit. In particular, the question as to what extent (if at all) such crossover 
scaling functions are universal is an intriguing one [9-11,22,23,33]. 

Another very interesting crossover which can also be studied is that which 
occurs near the critical point of unmixing for polymer solutions in a bad 
solvent [12,34-43]. For chain length ^ oo the critical temperature Tc{N) 
moves towards the G-temperature, where a single coil undergoes a transition 
from a swollen coil to a collapsed globule. This limit corresponds to a tricritical 
point [12]. 

Monte Carlo analyses of critical phenomena typically apply finite-size scaling 
concepts [44-49]. However, care is necessary in the proper application of these 
methods in the mean-field limit. In fact, the standard formulation of finite- 
size scaling ("linear dimensions L scale with the correlation length ^") im- 
plies that the hyperscaling relation [2] between critical exponents should hold 
[45,50], which is not the case for mean-field exponents (apart from d = d* = 4 
dimensions). This problem already arises for Ising models with short-range 
interactions for d > d* [51-61], and some disagreements between Monte Carlo 
results [51,52,54] and theoretical predictions [53,60] have stimulated a long- 
standing debate (see [61] for a detailed review). 



2 Mean-field to Ising crossover 

We consider the Hamiltonian [5,6] 

n/kBT = - E E ^i^i - - ^0 E ' (1) 

i j>i i 
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with Si = ±1 and an interaction K{r) = cR~'^ for |r| < R and zero elsewhere. 
The critical behavior of this model on rf- dimensional lattices can be studied ef- 
ficiently with a new cluster algorithm adapted for long-range interactions [62] . 

To analyze the crossover it is instructive to consider the associated Ginzburg- 
Landau field theory in continuous space, 



n{4>)/k^T^- 




(2) 



where 0(r) is the single-component order-parameter field, v is a temperature- 
like parameter and uq is a constant. After Fourier transformation and suitable 
rescaling this can be rewritten as (here N is the total number of lattice sites) 
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where u is related to Uq and h to Hq [6], and Tq in mean-field theory is the 
deviation of the temperature from its critical-point value. 

We are now interested in identifying the crossover scaling variable associated 
with the crossover from the Gaussian fixed point u = and ro = to the 
nontrivial Ising fixed point (Fig. 1). Because of the trivial character of the 
Gaussian fixed point and the fact the crossover scahng description should 
hold all the way from the Ising fixed point to the Gaussian fixed point, one 
can infer the crossover length scale Iq = R'^/(^~^^ exactly! This is done by 
considering a renormalization by a length scale /, such that the wavenumber 
changes from k to k' = kl, the number of degrees of freedom is reduced from 
N to N' — Nl~'^, and ■^^ changes into ■0^, = I'^i/^k to leave H invariant. From 
inspection of the terms in the Hamiltonian one can conclude that the singular 
part of the free energy must satisfy the scaling relation 

~ (ro_ u_ h\ _j-df ^j4-d Ki+dA 

^'[r'^'r^'r ~ ^'[r^ ' r^ ' r ■ ^ ^ 



We see that a finite and nonzero value for the second argument of fg is retained 
exactly when I takes the value of the crossover scale Iq. Thus we conclude that 
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the singular part of the free energy scales with R as follows 



(5) 



where a natural choice of coordinates (Fig. 1) is to measure fo and uq as 
distances from the Ising fixed point, unlike in the original Hamiltonian, where 
To and Uq are distances from the Gaussian fixed point. 

Equation (5) describes how the temperature distance fo from criticality and 
the magnetic field h scale with the range of interaction R: Note that the 
crossover exponent is known exactly (unlike other cases of crossover, e.g., be- 
tween the Ising and Hciscnbcrg universality class in isotropic magnets with 
varying uniaxial anisotropy [63]). The same result for the crossover exponent 
follows [5], of course, from simple-minded arguments using the Ginzburg cri- 
terion. However, the location of the nontrivial fixed point u* (Fig. 1), the 
associated other exponents, and the explicit form of the scaling function /g 
cannot be obtained exactly. 

The calculation of the scaling function for the free-energy density or its deriva- 
tives, such as the susceptibility, is a nontrivial task for both renormalization- 
group and Monte Carlo calculations. This is demonstrated in Fig. 2 where 
the effective critical exponent of the susceptibihty for T > Tc is plotted 
versus the thermal crossover scaling variable t/G, with t = (T — Tc)/Tc being 
the reduced temperature and G = GqR~^ the Ginzburg number in 0? = 3, for 
which Go ~ 0.277. Note that effective exponents are defined as 

7± = -d\nx/d\n |t| , X = WR){dM/dh)T , (6) 

where ± refers to T ^ T^, respectively, M = (s)t,/i) and the range R is defined 
from the second moment of the interaction [z being the effective coordination 
number) 



R' = Y: \r^ - r,\'K{v, - r,)/ ^ K{r, - r,) 

= - m - with |rj -Tjl <Rm- (7) 

Z ■ / ■ 

Here the second equality holds only for a square-well potential and values 
Rl^ = 1, 2, 3, 4, 5, 6, 8, 12, 18, 28, 60, 100, and 160 were studied. From Fig. 2 we 
see that the Monte Carlo results agree with all the theoretical calculations 
near the Gaussian fixed point, but do not yield the more rapid increase of 7^ 
near the Ising fixed point. It is not clear what precise conclusions should be 
drawn from this discrepancy: All these theoretical treatments really rely on 
extrapolations of low-order renormalization-group expansions in e — A — d, and 
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hence are perhaps rather inaccurate in d — 3 dimensions. On the other hand, 
they clearly relate to the limit where R ^ oo and t ^ 0. with tR^ fixed — a 
universal description of the crossover can only be expected in this limit. The 
Monte Carlo data shown in Fig. 2 also include the range of small R, for which 
additional corrections to scaling present near the Ising fixed point (other than 
those attributable to the Ising-mean-field crossover) may come into play. 

A rather successful description of the Monte Carlo data could be obtained by 
a fit to a function given by Anisimov et al. [33]. Their description is also an 
interpolation formula based on low-order e-expansions but contains a second 
parameter (in addition to G) describing a short-wavelength cutofi'. However, 
one disturbing feature of this description is that one needs different amplitudes 
Go in the relation G = GqR~^ above and below Tc, and the ratio Gq /Gq is an 
additional, ad hoc, parameter the significance of which is not understood [33]. 
Thus we consider it an as yet unsettled problem as to just on which parame- 
ters the crossover scaling description should depend. In this context, we draw 
attention to the question whether the specific square-well form chosen for the 
exchange interaction matters. To answer this question, a more general form 
of K{ri — Tj) was chosen (viz., a superposition of two square- well potentials 
which differ in range and strength but are chosen such that the same value 
for R^ results [10]. While was shown not to be determined by R alone, but 
depended on K{ri — Vj) in a more detailed way, the same crossover scaling 
function resulted for all choices of the interaction profile studied [10]. 

A particular merit of the description of Anisimov et al. [33] is, however, that 
it can yield a non-monotonic variation of 7<J^ with t/G: In d = 3 a shallow 
minimum (7^^ pa 0.96 < 7mf = 1) occurs for \t\R^ = 10^ [9] that can be 
fitted by this theory [33]. Indeed, a very similar minimum has been observed 
in Ref. [64] from a mean-field expansion for Ising systems with medium-range 
interactions, see also Ref. [65] for a detailed review. In c? = 2 dimensions, such 
a minimum occurs as well and is much more pronounced than in ci = 3, while 
above Tc the variation of the effective exponent is still monotonic (Fig. 3). 
Note that the crossover is again spread out over many decades in the crossover 
variable t/G {G oc in d = 2), as in d = 3, and that for T < there are no 
analytical results whatsoever to compare our Monte Carlo results with! At this 
point, there is clearly still a gap in our knowledge about critical phenomena. 

As a last point in this section, we add a few brief comments about the way in 
which the Monte Carlo results on effective exponents have been obtained. As is 
well known [45-49] , the Monte Carlo method converges to the exact statistical 
mechanics of a finite system only; the thermodynamic limit is never addressed 
directly. The typical situation is that one deals with slLxLotLxLxL 
box with periodic boundary conditions. The critical singularities are rounded 
and shifted by the finite size of the system [44-49]. For the precise location 
of the critical point, a finite-size scaling analysis is required. The principle of 
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finite-size scaling is that the hnear dimension L scales with the correlation 
length ^. Therefore the /c'th moment of the magnetization m scales like: 

{\M\^)^L-^^'''m{L/^), (8) 

(3 and v being the critical exponents of the order parameter ((m) oc \t\^) 
and the correlation length oc respectively, and Mk being some scal- 

ing function. Therefore the straightforward observation results [45] that these 
power law prefactors L~^^/^ cancel out if one considers suitable ratios of mo- 
ments, such as 

Q = {mYI{m') = g(L/e) . (9) 

At Tc we have ^ — > oo, of course, so (3(0) is simply a constant, indepen- 
dent of the system size L. This justifies the simple recipe to record this ratio 
for different choices of L and obtain from the intersection point of these 
ratios [45,48,49]. Note that the ordinate value of this intersection point is 
a universal constant (only depending on the shape of the system and on the 
boundary conditions, but not on i?, for instance, provided one is in the asymp- 
totic critical region). 

However, this recipe so far ignores the crossover from one universality class to 
the other (as well as corrections to scahng). Nevertheles, it turns out that one 
can formulate a combined finite-size scaling and crossover scaling description 
for such problems [5-11,16,17,66]. A simplified description considers the vari- 
ation of the correlation length, which is ^ oc Rt^^/'^ in the mean-field critical 
region, and ^ oc {R'^t)~'^ in the Ising critical region [the exponent k follows 
from the condition that for t = tcross oc i?~2<^/(^~'^) and the corresponding value 
of Across — ^{t — Across) = ^0 cx^ i?4/(4-(i) ^ smooth crossover between both 
power laws occurs]. Now it is of crucial importance to compare L with the 
crossover length scale Across- If L is much less than .^cross) then the finite size 
rounding occurs fully in the mean-field regime, before the crossover to Ising 
criticality has had a chance to come into play. Actually in this regime the 
correlation length ^ is not the relevant length to describe the finite size round- 
ing [48,51,52], one rather needs the so-callcd "thermodynamic length" [52], 
£t oc as will be discussed in Sec. 4. In this regime (L <^ Across) an 

accurate determination of Tc is clearly impossible. In order to accurately lo- 
cate Tc, we need to study the inverse regime, L 3> Across: Only then can one 
see the mean-field critical behavior farther away from Tc crossing over to the 
Ising behavior at tcioss (remember that this crossover is spread out over several 
decades!) and the finite-size rounding sets in at a still much smaller value of |t| 
(where L ~ ,^). Since for large R, .^cross is also very large (.^cmss — /o oc R'^^^'^'^^), 
one needs to simulate very large L and hence such simulations are technically 
very difficult. 
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Thus it is not surprising that when this problem was first addressed with 
single-spin-flip Monte Carlo algorithms [5] a satisfactory description of the 
full crossover could not be obtained, and the availability of an efficient cluster 
algorithm [62] was crucial for obtaining meaningful results. In d = 2, we could 
study L up to 800 lattice units, and Rm — 100 corresponding to z — 436 
interacting neighbors (Fig. 4). With these large lattices it is possible to follow 
the variation of Q almost all the way from the mean-field limit at small L to 
the Ising limit at large L, and in x (Fig- 4(a) the Ising asymptote (slope 3/4 
on the log-log plot) is nicely confirmed). 

Since we know T^. very precisely and have data for such a wide range of L, it is 
also possible to carry out runs slightly away from T^,, which are used to study 
the thermal crossover presented in Figs. 2, 3. Only data not affected by the 
finite system size are used for the numerical derivative required in Eq. (6). 



3 First steps towards the study of crossover problems in polymer 
blends and solutions. 

As is well known [3,12-14], the Ginzbiirg-Landau-Wilson Hamiltonian for a 
symmetrical polymer mixture near its critical unmixing point can be mapped 
on to the Ising model with a medium range of interaction (in o? = 3 dimen- 
sions), N^/'^ (with A'^ being the chain length of the flexible macromolecule) 
playing the role of the interaction volume B?. Qualitatively, this mapping is 
understood from the fact that a polymer coil has a random walk-like conflgu- 
ration. Its gyration radius Rgyr scales as Rgyr — N/Q^ where a is the size of 
the monomer. Thus the monomer density of one chain inside the volume that 
is occupies {V oc R^gyr) is very small, p = N/V oc a~^iV~^/^. Hence in a dense 
melt (pmeit — «~^) there are N^^'^ chains in the same volume, i.e. each chain 
interacts with x — N^^"^ "neighbors". Thus as — > oo one again expects a 
crossover from Ising-like critical behavior to mean-fleld like behavior, and this 
is verifled experimentally [18,19] (though the corresponding prediction for the 
Ginzburg number G (x 1/N does not seem to work out. 

First steps to study this crossover by computer simulation have been per- 
formed [16,17] using the bond fluctuation model of symmetrical polymer mix- 
tures [9,67] applying a semi-grand canonical algorithm [15] and histogram 
reweighting techniques [68]. The model and methodology of these simulations 
have been extensively reviewed elsewhere [3,67] and hence we omit all the 
technical details here, and simply show an attempt to estimate the crossover 
scaling function of the order parameter [17] (Fig. 5). Note that polymer are 
slowly relaxing objects and hence difficult to simulate — no counterpart to the 
cluster algorithm used for the Ising model [62] is available, and hence the chal- 
lenge remains to improve substantially the accuracy of studies such as shown 
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in Fig. 5 in order to be able to study the variation of effective exponents for 
this problem in analogy with Figs. 2, 3. 

If one wishes to compare such simulations for polymer mixtures to experi- 
ments on real systems [18,19], an important complication that must be taken 
into account is the asymmetry in chain length, Nj^ ^ Nb- This leads to two 
very important technical complications: (i) While in the symmetrical case the 
coexistence curve (including the critical point) occurs at a chemical potential 
difference A// = 0, for Na 7^ Nb phase coexistence occurs along a non trivial 
curve A/icocx(^) in the (A//, T) plane, and hence one has to search for the 
critical point (Aficrit = A/icocx(^c)i ^c) in a two dimensional variable space. 
Fig. 6. shows that this problem can also be overcome by finite-size scaling 
methods, utilizing the scaling behavior appropriate for first-order transitions 
AfjL— AncoexiT) oc L""*^ in d dimensions [48,49] in order to locate A//coex(r)[20]. 
(ii) Owing to the asymmetry, order parameter density and energy density be- 
come coupled, and this "field mixing" effect needs to be disentangled from 
the finite-size scaling analysis [21]. This problem is well known from computer 
simulation of fiuids and we shall not describe it here, but rather draw attention 
to a recent review [69]. 

This field-mixing problem is particularly severe for the unmixing of polymers 
in solution beneath the 0-temperature (which formally can be considered as 
a limiting case of a polymer mixture where Nb = N, Na = 1[3]), Fig. 7. How- 
ever, by a suitable transformation of variables, one can construct from and 
the energy density u, an appropriate field Ai = {(f) — su)/{l — sr) (where s,r 
are parameters that can be found from a suitable analysis of the simulations, 
see [21,69]), which then scales like the magnetization of the Ising ferromag- 
nets. Figure 8 shows that the distributions of this variable at criticahty nicely 
coincides with the critical order parameter distribution of the Ising model (ac- 
tually this mapping can be used as a method for precisely locating the critical 
point [21,69]). 

From analyses of this kind it has been possible to obtain the critical parameters 
of the model as a function of chain length, sec e.g. Fig. 8. The simulation 
results reproduce nicely the behavior oc A^^^' with x ^ 0.37 found also 
experimentally [41,42]. However, the simulations also show that the chains 
at the critical point are not yet partially collapsed, but are rather ideal, and 
hence rule out the interpretation of this exponent value (which differs from 
the classical results x — 1/2 [41]) as being due to the percolation of partially 
collapsed chains. Consequently, the physical interpretation of this exponent 
remains an open question [34,43]. 
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4 Finite-size scaling above the upper critical dimension 



Remembering that the correlation length ior d > d* — A has the mean-field 
critical behavior ^j, — ^ot~^^'^: the free-energy density can be written as [70] 

h-L-'fi^tl^^J ,uL^-^hL'^'/'^ . (10) 



Note that here exactly the same powers of L appear as those for / in Eq. (4) . For 
d > d* there is only the Gaussian fixed point to be considered. But although 
u* — here and the power of L in the term uL'^~'^ is negative so that uL^"'^ — > 
for L — > cxo, the argument uL^~'^ must not be omitted: m is a "dangerous 
irrelevant variable" [71], so when we consider the zero-field susceptibility x 
and the moment ratio Q [Eq. (9)], we find, using u oc £q~^, 

(12) 




Thus all scaling functions have two arguments, t{L/^o)'^ and {L/io)^~'^. How- 
ever, it turns out [51] that a reduction to one-variable scaling occurs for 
L — > oo, namely 

X - Jim L^/^P, (tL'^/%^£t'^^') , (13) 
Q ^ hm Pg {tino'^t'^^") = pQiiL/err/'} , (14) 



where in the last step we have introduced the "thermodynamic length" i-r oc 
]^^2], mentioned above. 

Equations (13,14) can be understood from various arguments [51-53]. Brezin 
and Zinn- Justin argue [53] that in the initial Hamiltonian or the corresponding 
statistical weight, one can treat the contribution from the average magnetiza- 
tion M separately, 

«xp[-«{..}ABrl^«xp{-(M^|^L^ + ...} . (15) 



where the dots stand for contributions with non-uniform magnetization, i.e. 
fiuctuations. Here Mb, Xb are the mean-field bulk magnetization and suscep- 
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tibility, Mh = Mb{—ty^'^,Xb = Xftl^l"^- The zero-mode theory neglects these 
fluctuations altogether and there the distribution of the magnetization Pl{M) 
scales as 



Prom this result it is straightforward to derive the above scaling functions 
and pQ explicitly [53]. 

Since this theory was proposed [51-53] it has been a long-standing problem 
to verify the predictions by Monte Carlo simulation. In particular, when one 
plots the moment ratio Q versus temperature deviation from criticality, one 
should find a universal intersection point at T^. at a value 



However, the Monte Carlo results for small systems seem to intersect at a 
different value Q ~ 0.52 (Fig. 10). Also the temperature where this intersection 
occurs is a httle off, but since one does not know Tc beforehand, one could 
simply imagine that the abscissa in Fig. 10 is mislabeled and must be 
assigned differently. 

Chen and Dohm [59,60] have recently criticized the whole approach sketched 
above and maintained that one must return to a finite-size scaling descrip- 
tion in which both variables t{L/^)^ and {L/IqY~'^ are kept separate, as in 
Eqs. (11) and (12). They also obtained the scaling functions Pq and in a 
first-order loop expansion as a function of these variables. Indeed their result 
is quahtatively similar to the Monte Carlo data (Fig. 10, broken curves), al- 
though in quantitative respects their treatment offers little improvement. This 
is seen, for instance, in a scaling plot of the susceptibility: The Chen-Dohm 
theory approaches the zero-mode results from above, while in the regime of 
interest the Monte Carlo data fall below the zero-mode result (Fig. 11). These 
discrepancies remain present for considerably larger L than shown here [61]. 

Thus we arrive at a rather disappointing state of affairs — although for the 0^ 
theory in o? = 5 dimensions all exponents are known, including those of the 
corrections to scaling, and in principle very complete analytical calculations 
are possible, the existing theories clearly are not so good. Perhaps the discrep- 
ancies result because the theory of Rcf. [60] is only one-loop order, perhaps 
because other corrections arc missing. While presumably the zero-mode one- 
parameter scaling is true asymptotically for L — oo, the corrections to this 
limit disappear only rather slowly, as Fig. 10(a) has demonstrated. 



Pl{M) oc L''l^eM-[M^I{Ml{-t)) - IWL/ItYI^} ■ 



(16) 



Pq(0) = 87rVr^(l/4) ~ 0.456947 . 



(17) 
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5 Concluding remcirks 



While the estimates of the critical exponents for the d ~ 3 Ising model are 
impressively accurate [72-74] and analytical [72] and Monte Carlo [73,74] esti- 
mates agree within very small error margins, the situation is different for the 
problems considered in the present paper: Analytical work is restricted to low- 
order ^-expansions or low-order loop-expansions and discrepancies between 
theory and simulation occur that are not fully understood. More work will be 
needed to clarify the situation. Note that the Ising to mean-field crossover con- 
sidered here really is the simplest example of crossover phenomena, since the 
crossover exponent is rigorously known — crossover from one nontrivial fixed 
point to another is presumably more tricky to deal with. And for problems 
such as the critical point of polymer solutions, even the proper theoretical ap- 
proach is controversial, and hence it is unclear whether the exponent x ~ 0.37 
(Fig. 9) is a universal property at all [34-43]. 

Further problems appear when one is not concerned with bulk critical phenom- 
ena in ideal, homogeneous systems, but when one considers inhomogeneous 
systems, e.g., systems with random quenched disorder (e.g., Ising and Potts 
models exposed to random fields, spin glasses, etc. [75]). For instance, for a 
Potts spin glass finite-size scaling is not even understood on the mean-field 
level, at least for cases where first-order transitions without latent heat oc- 
cur [76]. Also for systems with a regular inhomogeneity, e.g., Ising films with 
competing walls which allow for interface localization-delocalization transi- 
tions, one has fascinating critical behavior and crossover, of which the details 
still need to be unraveled [77] . Thus the Monte Carlo investigation of phase 
transitions — both in equilibrium and in driven systems [78] — will remain an 
active and challenging field. 
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Fig. 1. Qualitative picture of the renormalization trajectory describing the crossover 
from the Gaussian fixed point //q = (0,0) to the Ising fixed point fi* = {rQ,u*). 
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Fig. 2. Effective susceptibility exponent 7*^ above Tc for the three-dimensional 
Ising model with variable interaction range R (numbers in the key) plotted vs. t/G, 
along with three theoretical calculations for this quantity; due to Refs. [31] (BK), 
[27] (BB), and [25,29] (SF), respectively. From Ref. [9]. 
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Fig. 3. The effective susceptibility exponent 7*g above Tc (a) and below Tc (b), for 
the two-dimensional Ising model with variable interaction range R (numbers in the 
key), plotted vs. tE?. Prom Ref. [8]. 
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Fig. 4. Finite-size crossover curve for the magnetic susceptibility x divided by 
the system hnear dimension (a) and the amplitude ratio Q [Eq. (9)] (b) for the 
two-dimensional Ising model a.t K = KdR) plotted vs. the finite-size crossover 
scaling variable L/B^ (note that Across = E? in d = 2). In both quantities, 

range-dependent correction factors C[x] and C[Q] have been divided out to elimi- 
nate some corrections to scaling (see Ref. [7] for a definition of these factors). Prom 
Ref. [7]. 
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Fig. 5. Crossover scaling plot for the order parameter (|m|) = {\4)a — (t>B\) / [4>a + 4>b) 
of a binary polymer mixture (A,B) with symmetrical chain lengths A^^ = Nb = N . 
(l>A,(pB are the volume fractions of A and B monomers, respectively. The points are 
simulation results for the bond-fluctuation model on a simple-cubic lattice, using 
concentration cpy = 0.5 of vacant sites. Straight lines in this log-log plot indicate 
power laws with effective exponents, (m) = B^st^'^^ , t = 1 — T/Tc- The broken 
straight line shows the mean-field result, (m) = y/St^^"^, to which the data converge 
for N ^oo. Prom Ref. [17]. 
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Fig. 6. Finite-size scaling plot for the second moment of the order parameter of 
an asymmetric polymer mixture (iV^ = 10, Nb = 20) at a temperature T < 
{e = £ab/T = 0.035) as a function of the normalized chemical potential difference, 
in order to locate A/Xcoex(T') by optimizing the "data collapse" for the range of 
values of L as indicated. From Ref. [20]. 
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Fig. 7. Schematic phase diagram of a polymer solution using the temperature T and 
the volume fraction taken by the effective monomers of the polymer chains as vari- 
ables. The coexistence curve separates a dilute solution of collapsed chains (at i^iolx) 
from a scmi-dilutc solution of overlapping chains (at 0cocx)- These two branches of 
the coexistence curve merge at a critical point Tc{N), (pc(N). For N ^ oo this point 
merges with the G-point of a polymer solution at infinite dilution — ^ 0). 
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Fig. 8. Critical order-parameter distribution for a polymer solution with chain length 
N = 20, modeled by the bond-fluctuation model on the simple-cubic lattice, for 
linear dimensions L = 40 and L = 50, open symbols, and compared to the or- 
der-parameter distribution of the Ising model (crosses). Prom Ref. [43]. 
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Fig. 9. Estimates of the critical volume fraction (pc of monomers for a poly- 
mer solution (modeled by the bond-fluctuation model on the simple-cubic lattice, 
with an attractive energy between monomers at distances r < \/6) as a func- 
tion of the inverse chain length. The broken curven represents a fit of the form 
= (1.1126 + 1.3A^°-369)-i. Prom Ref. [43]. 
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Fig. 10. (a) Plot of Q vs. tL^ for the d = 5 Ising model, demonstrating the occurrence 
of spurious intersections both in the Monte Carlo results [61] and the Chen-Dohm 
theory [60]. (b) Plot of Q vs. the scaling variable tL^/"^; using parameters 
extracted from various limits of the susceptibility [61], the Chen-Dohm theory can 
be evaluated without any adjustable parameter whatsoever. Note that for L = 12 it 
is already graphically indistinguishable from the "zero-mode" theory. Prom Ref. [61]. 
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Fig. 11. Plot of the scaled susceptibility vs. tL^^'^, including the "zero-mode" 

result of Ref. [53] , as well as the predictions of Chen and Dohm [60] evaluated for 
the same values of L as the Monte Carlo results shown. Prom Ref. [61]. 
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